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Abstract 

It is well known that steady states of the Vlasov-Poisson system, a widely used model in 
non-relativistic galactic dynamics, have negative energy. In this paper we derive the analogous 
property for two relativistic generalizations of the Vlasov-Poisson system: The Nordstrom- 
Vlasov system and the Einstein- Vlasov system. In the first case we show that the energy of 
steady states is bounded by their total rest mass; in the second case, where we also assume 
spherical symmetry, we prove an inequality which involves not only the energy and the rest 
mass, but also the central redshift. In both cases the proof makes use of integral inequalities 
satisfied by time depedent solutions and which are derived using the vector fields multipliers 
method. 
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1 Introduction and main results 

A widely used model in astrophysics for the dynamics of the stars of a galaxy is the Vlasov-Poisson 
system [7]. This model is justified when collisions among the stars and external forces are neglected. 
Moreover, being a non-relativistic model, the Vlasov-Poisson system ceases to be valid when the 
stars move with large velocities (of the order of the speed of light) or in the presence of very 
massive galaxies, since then relativistic effects become important. Typical relativistic effects are 
the redshift of the luminous signals emitted by the galaxy and the formation of black holes. The 
model which is currcntely believed to represent the physically correct relativistic generalization of 
the Vlasov-Poisson system is the Einstein- Vlasov system [T] , where Poisson's equation is substituted 
by Einstein's equations of General Relativity. Another relativistic generalization of Vlasov-Poisson 
is the Nordstrom- Vlasov system [5]. This model relies on the same geometric interpretation of 
gravity as for the Einstein- Vlasov system. Although it is not physically correct, the Nordstrom- 
Vlasov system is mathematically interesting since it already captures some of the technical and 
conceptual new difficulties that are encountered when studying a relativistic (Lorentz invariant) 
system. 

In this paper we want to investigate the mass-energy bounds (virial inequalities) required for the 
existence of steady states to the relativistic models. It is well known in fact that static solutions 
of the Vlasov-Poisson system, which correspond to equilibrium configurations of the galaxy, have 
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negative energy. The same property cannot of course be true for the relativistic models, since the 
energy in the latter case is always positive. Before considering the relativistic models in detail we 
present the role of the virial identities in the case of steady states to the Vlasov-Poisson system. 



1.1 The classical case: Vlasov-Poisson system 

Let / = f(t, x,p) be the distribution function in phase space for an ensemble of unit mass particles, 
where / > 0, t G K, x € K 3 and pgl 3 . In the physics of gravitational systems, the particles stand 
for the stars of a galaxy. In geometric units, i.e., AttG — 1, where G is Newton's gravitational con- 
stant, the gravitational potential U — U(t,x) generated by the galaxy solves the Poisson equation 

A x U = p, lim U = Q,Vt€R, (la) 

| IE | — >00 

where 

P= I fdp (lb) 



is the mass density of the galaxy and the boundary condition at infinity means that the galaxy is 
isolated. The assumption that the stars interact only by gravity leads to the Vlasov equation: 

9*/+p-V x /-V x C/-V p / = 0. (lc) 

The system ([T]) is the Vlasov-Poisson system. The energy H and the mass M of a solution are 
given by 

H = ~ [ I \p\ 2 fd P dx-\ f \V x U\ 2 dx, M= [ f fdpdx 
1 JR 3 JR 3 1 Jr 3 Jr 3 Jr 3 

and are conserved quantities. Likewise, the total linear momentum Q and angular momentum L, 

p f dpdx , L = / / x Ap f dpdx , 



are conserved quantities. The invariance of Vlasov-Poisson by (time dependent) Galilean transfor- 
mations is the property that, given !i£l 3 and the transformation of coordinates 

Q u : t' — t , x' — x — ut , p' — p — u , 

then f u (t, x,p) — f(t' 7 x\p') and U u (t, x) = U(t', x') solve the system (J]) if and only if (/, U) does. 
Note that Q can be made to vanish by a Galilean transformation with velocity u — Q/M; the 
resulting reference frame is at rest with respect to the center of mass of the distribution, which is 
defined as c p (t) = M" 1 J R3 x pdx . 

A galaxy in equilibrium is described by steady states solutions of the Vlasov-Poisson system. We 
distinguish between two types of steady states: Static solutions and traveling steady states. The 
formers are defined as time independent solutions of the Vlasov-Poisson system JT]) and have total 
momentum Q — 0. A solution / is a traveling steady state (with total momentum Q =/= 0) if 
/ ° Gut where u — Q/M, is a time independent solution of the Vlasov-Poisson system (i.e., a 
static solution). Our interest on traveling steady states is motivated by the fact that their energy 
provides a lower limit for the energy of totally dispersive solutions, see Moreover, the non- 
linear stability theorems proved for the Vlasov-Poisson system consider the traveling steady states 
as possible perturbations of a static equilibria, see [TOj, [19l [24] and references therein. 



A fundamental property shared by all static solutions of the Vlasov-Poisson system is that of having 
negative energy. The proof goes as follows. Any sufficiently regular solution of the Vlasov-Poisson 
system satisfies the dilation identity: 

-77 / / x ■ p f dpdx = H + E kin , E kin = i / / \p\ 2 fdpdx, 

at JR 3 JR 3 Z JR 3 JR 3 
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as it follows by direct computation. If / is a static solution, then the previous identity implies the 
virial relation H = — £"kin, which yields that 

H < , for static solutions of the Vlasov-Poisson system. (2) 

For traveling steady states, we just apply to {2J a Galilean transformation with u = (M)~ 1 Q and 
we obtain 

Q|2 

H < i f° r traveling steady states of the Vlasov-Poisson system. (3) 

Our purpose in this work is to extend these fundamental inequalities to the relativistic case. For 
more information on the Vlasov-Poisson system, we refer to [HJ HH [17l Qjj] . 



1.2 Main results for relativistic models 

In the case of the Nordstrom- Vlasov system, which will be presented in Section[2l the generalization 
of ((21) is that the energy of regular steady states is bounded by their mass, i.e. 

H < M . (4) 

Furthermore, the counterpart to ([3]) for traveling steady states is 



(Of course, the invariants H, M , Q have to be redefined in an appropriate way for this new system, 
see Section [2]). Moreover the equality sign could only hold for steady states with unbounded 
support. We remark that in the case of the Vlasov-Poisson system the supremum of the steady 
states energy coincide with the infimum energy of totally dispersive time dependent solutions, 
see [TT]. The analogous statement for the Nordstrom- Vlasov system is currently not known, due 
to the difficulties in defining a Lorentz invariant concept of total dispersion. We also remark that 
the bound H < M, which holds for all regular and compactly supported static solutions of the 
Nordstrom- Vlasov system, is crucial in the proof of orbit stability of the polytropic steady states 
established in [ID] . 

For the spherically symmetric Einstein- Vlasov system, analyzed in Section (3[ we derive an inequal- 
ity that involves not only the energy (ADM mass, H) and the mass (rest mass, M) of the steady 
state, but also the central redshift Z c : 



Z r > 



M 

1 

H 



(5) 



The metric of the space-time for spherically symmetric static solutions of the Einstein- Vlasov 
system is determined, following the notation in Section [3l by two functions A(r) > and /i(r) < 
of the radial variable. The central redshift is defined in terms of the second one at the origin 
by Z c := e~''(°) — 1. It is the redshift of a photon emitted from the center of the galaxy. The 
estimate ([SJ can thus be seen as an upper bound for fj.(0). Similarly, the celebrated Buchdahl's 
inequality in General Relativity [27] can be seen as an upper bound on the metric component A(r) 
for spherically symmetric steady states of the Einstein-matter equations. A quite general version 
of the Buchdahl inequality was proved recently in [2] and reads 

/ \ 8 

sup (l — e _2A ^ r M < - , or equivalently supA(r) < In 3 . (6) 

For static shells the Buchdahl inequality is equivalent to a lower bound for the external radius. 
We will show that estimate (JSJ leads to an upper bound on the internal radius. We refer to [5] for 



3 



an analitical/numerical investigation of the Buchdahl inequality in the context of the spherically 
symmetric Einstein- Vlasov system. 

We do not know whether, as for the Vlasov-Poisson and the Nordstrom- Vlasov system, the inequal- 
ity © could also be related to the problem of stability of spherically symmetric static solutions. 
This is a difficult question to answer, since the stability problem for the Einstein- Vlasov sys- 
tem is still poorly understood. However it is worth noticing that heuristic and numerical studies 
[29l [30l [26] indicate that the regime of stability of compact galaxies is indeed characterized by 
the central redshift and the fractional binding energy (defined as 1 — H/M). Moreover it was 
conjectured that the binding energy maximum along a steady state sequence signals the onset 
of instability. There are several numerical studies on the problem of stability for the spherically 
Einstein- Vlasov system; we refer to [U [5] [55] . 

A last basic comment on ([5]) is that, as opposed to the inequalities that hold for steady states 
of the Vlasov-Poisson and the Nordstrom- Vlasov system, the bound §5§ contains a quantity, the 
central redshift, which is not preserved along time dependent solutions. It is therefore not clear 
whether one can interpret ([5]) as the exact analog of the mass-energy inequalities for the steady 
states of the Vlasov-Poisson and Nordstrom- Vlasov system. 

We conclude this Introduction with a brief explanation on how we prove our main results. As 
a first step we employ the vector fields multipliers method to the local conservation laws for the 
Nordstrom- Vlasov and the Einstein- Vlasov system to establish a virial identity which has to be 
satisfied by all time dependent solutions. These identities arc of independent interest and could 
be useful to derive space-time (Morawetz type) estimates for the evolution problem. The virial 
identities restricted to time independent solutions give rise, after applying some simple bounds on 
the moments of the distribution /, to the virial inequalities ([l])-©. 



2 The Nordstrom- Vlasov case 

We write the Nordstrom- Vlasov system in the formulation used in [TO]: 



ft/ + ^ |p|2 • V,/ - V, (y«* + W ) • V P / = , (7a) 

flfr-A.^-e*/ /_*_. (7b) 
JR3 y/eW + \p\ 2 

Here / = f(t,x,p) > and <j> — <j>(t,x). The physical interpretation of a solution (/,</>) is the 
following: The space-time is the Lorentzian manifold (R. 4 ,g = e 2 ^7y), where ?y is the Minkowski 
metric, whereas / is the kinetic distribution function of particles (the stars of a galaxy) moving 
along the geodesic curves of the metric g. The motion along geodesies reflects the condition that 
gravity is the only interaction among the particles. The system has been written in units such 
that AttG = c = 1, where G is Newton's gravitational constant and c the speed of light. If the 
latter is restored in the equations, in the limit c — -» oo one recovers the Vlasov-Poisson system 
in the gravitational case. For a proof of the latter statement and general information on the 
Nordstrom- Vlasov system we refer to [12] . The global regularity of solutions is studied in [8] [13] . 

The local energy, momentum and stress tensor of a solution (/, <f>) of JJJ) are defined respectively 
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as = 1,2,3) 

h(t,x)= f Je^ + \ P \2fdp+ l -{d t <t>) 2 + \\V^\ 2 
qi(t,x)= / pifdp — dt<i>di<j>, 



r ij (t,x) = [ |ff / dp + fti^ + \Sa [{M? - |V^| 2 ] , 

where di denotes the partial derivative along x l . These quantities are related by the conservation 
laws 

9th + V x ■ q = , 8 t q t + djTij = , (9) 

the sum over repeated indexes being understood. Upon integration, the previous identities lead to 
the conservation of the total energy and of the total momentum: 



H(t) = / h(t, x) dx — constant , Q(t) = / q(t, x) dx = constant . 

Moreover, solutions of the Nordstrom- Vlasov system satisfy the conservation of the total rest mass: 

M(t) = / p(t, x) dx = constant , 

which is obtained by integrating the local rest mass conservation law 

9 tP + V x -j = 0, p= [ fdp, j=l P fdp. (10) 
Jr 3 Jr 3 V« * + \P\ 

The system (J7]) satisfies the fundamental property of Lorentz invariance. Precisely, let (t',x') be a 
system of coordinates in Minkowski space obtained from (t, x) by a Lorentz boost, that is 

t = UqI — u ■ X , X = x — lit H ; — — [u ■ x)u , 

\W 

where u is a fixed vector in M 3 and uq = yl-j- \u\ 2 . The inverse Lorentz transformation is obtained 
by exchanging u with —u, that is 

t = uot' + u • x 1 , x = x' + ut' H — j (u ■ x')u , (11) 

\W 

which we shorten by (t, x) — L u (t', x'). Define the field <j) u in the new coordinates by 

<j> u (t',x') = 4>oL u (t',x'). 
Introduce the new momentum variable 

p = p — U\l 'e 2 ^*'*) + \p\ 2 H — " (u ■ p)u 

v ' \uY 

or, inverting, 

p=p' + uJeW*,*) + Ip/p + ( u ■ p')u . (12) 

v \u\ 2 

We shall write (t,x,p) = C u (t',x',p') to shorten the set of transformations (jlip - (fT2^) . Finally, 
define the distribution function in the new variables as 

f u (t',x',p') = fa£ u (t',x',p'). 
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In the language of special relativity, one says that / and <fi transform like scalar functions under 
Lorentz transformations. The Lorentz invariance of the Nordstrom- Vlasov system means that the 
pair (/,</>) solves the system ([7]) in the coordinates (t,x,p) if and only if (/«,</>„) satisfies the 
same system in the coordinates (t', x' ,p'). Thus, in particular, also the mass-energy-momentum of 
{fu, <Au) is conserved along the time evolution, 

M[f u ] = constant , H[f u , <p u ] — constant , Q[f u , <j>u] = constant . 

We shall need the relation between the mass-energy- momentum of (/, cf>) and of (/„, 4> u ), which is 
derived the following lemm£0. 

Lemma 1. For all u £ M 3 , 

M[f u ] = M[f] , (13a) 



H[f u , <p u ] = ^l + \u\ 2 H[f, 4>] - Q[f, 4>] ■ u , (13b) 

Q[U, <pu] = Q[f, 4>] - H[f, tf>]u + ^^(u ■ Q{f, <f>])u . (13c) 

\u\ 

Proof. Since the mass-energy-momentum of both pairs (/,</>) and (f u ,</>u) is conserved, it is suf- 
ficient to prove the relations (fl3|) for the initial value of M(u) :— M[f u ], H{u) := H[f u ,(j> u ] and 
Q(u) :— Q[f u , <Pu]- We restrict ourselves to prove the invariance of the total mass, the proof for 
the other transformations being similar. We shall need that, by (TH 



^J e 2<k(t,x) + | p /|2 = UQ J e 2<t>(t,x) + | p |2 _ u . p 

or, inverting, 

J#WM) + \p\ 2 = uo^je 2 ^) + \p'\ 2 +u-p. (14) 
In order to prove (|13a|) we write 



M{u)= / f u (0,x',p')dx'dp' = / foC u (0,x',p')dx'dp' 
Jr 3 Jr 3 Jr 3 Jr 3 

= 11 f(u-x',x' + ^-^(u ■ x')u,p' + U «/ e 2M0,z') + |p'|2 + r ^JZl( u . j/) u \ dx ' dp'. 

Jr 3 Jr 3 V M v M 2 J 

Next we make the change of variable 

x = x' H — 7 „ (u- x')u , p = p' + u\l e 2 ^"^' 1 ') + Ip'l 2 H — (u ■ p')u . 
\uy v \u\ z 

The Jacobian of this transformation is given by 



_u-p' + v / e*M + \p'\ 2 

v / e 20 u (O, a; ') + | p /|2 

where u = u/uq. Using (fl~2j) and (fT4)) we obtain 



J= 1- U ' P 



1 In the language of special relativity, the lemma establishes that M transforms like a scalar function, whereas 
the quadruple (H, Q) transforms like a four-vector under Lorentz transformations. 
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Since the volume measure transforms as dx'dp' = J 1 dxdp, we obtain 

M(u) = I I f(u-x,x,p) I 1 = ] dxdp 

JR 3 JR 3 \ , / e 24>(u-x,x) + \ p \2 

(p(u ■ x, x) — u ■ j(u ■ x, x)) dx , 

where p and j are defined by ()10[> . Taking the partial derivative d Ui of the previous expression we 
get 



d Ui M(u)= j (d t pd Ut (u ■ x) - (d Ui u k )j k ~u k d t jkd Ui (u ■ x)) (u ■ x,x)dx 
(-d Ui (u ■ x){d x J k + u k dtjk) - (d Ui u k )jk) (u ■ x, x) dx 
(d Ut (u ■ x)d Xk [j k (u • x, x)} - {d Ui u k )jk{u ■ x, x)) dx 

JR 3 

where we used the continuity equation (|10p to pass from the first to the second line and integration 
by parts to pass from the third to the last line. Thus we obtained that V u M(u) — 0, i.e., 
M(u) — M(0), which yields the claim on the invariance of the total rest mass. □ 

Remark 1. According to the transformation law of the total momentum Q, the Lorentz transfor- 
mation that makes Q to vanish, i.e. that moves the reference frame to the center of mass system^), 
is the transformation C u with u = Q/ \]H 2 — |Q| 2 . The energy of the transformed solution is 

VH 2 -\Q\ 2 - 



2.1 Virial identities for time dependent solutions 

The conservation laws can be expressed in a more coincise form as 

du,T1 = 0, n,u = 0,...,4, x°=t, (15) 

where is the stress-energy tensor, whose components are given by 

Too = — h , Tqi = —qi , T. t j = Ty . 

Indexes are raised and lowered with Minkowski's metric 77 M „ = diag(— 1, 1, 1, 1). Upon multiplying 
the conservation law (fT5j) by a vector field £ M = x), integrating on a compact spacetime region 
f2 with piecewise differentiable boundary dQ and applying the divergence theorem we obtain the 
integral identity 

/ TtOv^= [ T^d^dtdx, (16) 
Jan Jn 

where rt M denotes the exterior normal vector field to the boundary <9f2 and da the invariant volume 
measure thereon. The identities obtained from (fTrj|) upon a specific choice of the vector field 
multiplier go under the general name of virial identities. We prove here one that applies to regular 
asymptotically flat solutions. By this we mean that / 6 C , 4> € C 2 l~l L 2 , the mass and energy are 
finite and 

lim / h(t, x)dS R , = , VtER. (17) 

-R^oo Js(R) 

2 To be more precise, this is called the center of momentum system. In relativity there is no general acceptance 
on the concept of center of mass. 
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By oj we shall denote the outward unit normal to S(R) — [x : \x\ — R}, and (ISr stands for the 
invariant volume measure on S(R). Moreover we denote by x( r )i r > 0, a function that satisfies: 



x e C 2 , x' e L°° , — £ C 2 n L°° 



(18) 



Lemma 2. Let 



Z(t)= / x{r)(q-u-r 1 <j>d t cj>)dx , r = \x\ . 
For all regular asymptotically flat solutions of ([7]) the following identity holds: 



dl 



— = / x hdx+ -e 
at Us ins r 



X 2<f>i 



f 



: dpdx 



V r 



X 



\u A ~S7 x (j)\ + 



R3 ^e 2 ^ + \p\'- 



, \u Ap\ 2 + e 2 ^ 



y/e 2 <t> + |p| 2 



fdp 



dx 



^ —6 2 dx . 



2 Jrs r 

Proof. In (m]) we use Q = [0,T] x B(R), where = {x : \x\ < R} and 

e ■■ e = o,e = x(»v • 

We obtain 



(19) 



x(r)<?-wda; 



B(-R) 



X(r)TijU l u j dS R dt 



o JS(R) 
T 



Jfl(fl) 



dx dt 



(20) 



where for any function g(t) we denote [^(t)]^ = g{T) — .9(0). Using the bound \TijLJ l uP\ < 3h 
and lfT7|) we get 



X^TywVdSfl 
Then, letting i? — > oo in (j2"0)) we obtain 



<3||x||oo / hdS R ^0, R -> oo 

's(Ji) 



o -/o 



rJ r 



whence 



We compute 



and 



dt 



X(r)q -ujdx 



Fry = 



(x' - J) r yW V + 2f 



da; dt 



da; . 



8 



Hence 



4 / X{r)q -udx 
at j M 3 



3 2 2 



dx 



dx . 



Using that \ui A y\ 2 = |y| 2 — \uj ■ y\ 2 , for all vectors y € M 3 , we can rewrite the previous equation as 



^/ X (r)q-ojdx=f ^((d^-lV^dx- f X ' 



X 

— X 

r 



e^fdp 
w Ap\ 2 f dp 



dx 



Ve 2 * + |p| 2 

Moreover, using (|7bj) and integrating by parts twice, we find 



+ |wa v^r (21) 



4/ ^<pd t cpdx=( * ( (a t 0) 2 - |v^| 2 - ^ 



e 2 ^/ dp 
^/e 2 ^ + |p| 2 



A 



B(R) 



V r 



' dx + / —(j)uj-V x 4>dSji 

ls(R) r 



u) ■ V 



S(R) 



(7) 



(22) 



Applying the Cauchy-Schwartz inequality, the regularity of the solution and the assumptions on 
X, we obtain 



— <j>u ■ V x 4>dS R 
S(R) r 



< C\\(f>\\ L 2 {S{R)) \\V x (t)\\ L 2 {S{R)) 



<CJ h{t,x)dS R -► , R -> 00 , 



and 



w • V 



S(R) 



1? 



S(fl) 



C 
R 



where C is a constant independent from R. Thus taking the limit R — > 00 in (|2"2")) we get 



)d t (bdx 



Z ( (d t cp) 2 - |v^| 2 - ^ / ; e20/dp 1 <& 

I J R 3 ,/e 2 * + |p| 2 / 



The quantity A (^) is nothing but — . The sum of (|2~lj) and (|23[) yields the desired result. 



(23) 
□ 
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2.2 Virial inequalities for steady states 



As in the Vlasov-Poisson case, we distinguish between two types of steady states. Static solutions, 
which are defined as time independent solutions of the Nordstrom- Vlasov system ([?]) ■ and traveling 
steady states, which are defined as solutions f(t, x,p) such that / o C u , where u = Q/^/H 2 — \Q\ 2 , 
is a time independent solution of the Nordstrom- Vlasov system (i.e., a static solution). For static 
solutions one has Q = 0, whereas Q ^ for traveling steady states. Note that for static solutions 
(that vanish at infinity) the field is determined by / through a non-linear Poisson equation. Thus 
when we refer to a steady state solution we mean simply the distribution function /. The main 
goal of this section is to prove the following property of steady states to the system ([7]). 

Theorem 1. Let f be a static regular asymptotically flat solution of [[7]). Then 

H <M , (24) 

Traveling steady states satisfy y H 2 — \Q\ 2 < M. Moreover, equality in (|24p implies that the 
support of the static solution is unbounded. 



Proof. The statement on traveling steady states follows by applying the Lorentz transformation 
C u with u = Ql 'yj H 2 — \Q\ 2 to the inequality for static solutions, thus it suffices to prove the 
latter. To this purpose consider a function \ that, in addition to (|18[) . satisfies 



--X'>0, x"<0. (25) 
r 

Next we observe the simple inequality y — 1 > — e~ y , with equality if and only if y — 0. Using 

( f>-l>- e -' p (26) 

in the identity (Til?)) we obtain 

-?')*■ (27) 

In particular, for time independent solutions we have 

/ (x'h --p)dx<0. (28) 
Jr3 V r ' 

Let R > and consider the function x(r) = xr(t) given by 

( r for r ^ R , 

X{r) = \ 3R - 3^ + 4 for r > R . 

This function satisfies the properties ffl8]) and (|25|) . The left hand side of ([28|l becomes 

/ (x'h-^ P )dx= f (h-p)dx+f (x'h-^p) 
Jr^ r ' Jb(K) Jb{ry v r ' 



> (h- p)dx-C / hdx+ pdx] 

IB(R) \Jb(R)" JB{R) c J 

(h -p)dx + e{R) , (29) 

B(R) 

where s(R) — > as R — > oo. Thus, assuming > M, there exists i?o > such that e(i?) < 
(# - M)/4 and J B(R) (h - p) dx > (H - M)/2, for all i? > i? , whence 



/ (x'^--p)^>t(^-M)>0 
Jr3 V r / 4 
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which contradicts (|28|) . This concludes the proof of (|24|) . To prove that the strict inequality holds 
for static solutions with compact support, we observe that in the latter case the field <p never 
vanishes in the support of / (it is strictly negative) and thus the stronger inequality <j) — 1 > 
holds instead of ([26]) . Thus also the inequality in (|28]l is strict. Since the last member of ([29]) goes 
to zero for R — > oo when H = M, the claim follows. □ 

Remark 2. Theorem [T] improves a similar result proved in [TU] in two aspects. Firstly, in [TU] the 
fact that the strict inequality holds for compactly supported steady states was overlooked. Secondly 
the result presented here requires less decay than the inequality proved in [10] and therefore applies 
to more general steady states. In particular, this result allows to remove some technical hypothesis 
in the stability result obtained in [10] . 



3 The Einstein-Vlasov case 

The spherically symmetric Einstein-Vlasov system in Schwarzschild coordinates is given by the 
following set of equations (in units G = c = 1): 

d t f + e»- x jJ== • V,/ - (A t — + e^Vv'lThf ) - • V / = , (30) 
V 1 + M 2 V r It 



e~ 2X (2r\ r - 1) + 1 = 87rr 2 h , (31a) 

e- 2A (2r A t r + l)-l = 87rrV ad , (31b) 

A t = -47rre A+M g , (31c) 

' ~ 2A ( Mrr + (Pr ~ K){Hr + ~ e" 2 " (Aft + A t (A f - fl t )) = iltp^ , (31d) 



where 



h(t,r)= ^TTWfdv, p lad (t,r)= (—) f 

. B3 _ 1B3 \ T I 



, 2 „ dv 



x ■ v 



q{t,r) = fdv, p tan (t,r) 



r 

x A v 2 



dv 



The functions p rad and p tan are the radial and tangential pressure; h is the energy density and q 
the local momentum densit^. As usual, / > is the distribution function of particles (stars) in 
the phase space in the coordinates t € K, x £ M 3 , v £ R 3 . The variable v is not the canonical 
momentum of the particles, the latter being denoted by p in the previous sections. The function / 
is spherically symmetric in the sense that /(t, x, v) — f(t, Ax, Av), for all A £ SO(3). The symbol 
A denotes the standard vector product in R 3 , and for a function g = g(t,r), r = \x\, we denote 
by gt and g r the time and radial derivative, respectively. By abuse of notation, g(t,r) = g(t,x) 
for any spherically symmetric function. The functions A, p, determine the metric of the space-time 
according to 

ds 2 = ~e 2 ^dt 2 + e 2X dr 2 + r 2 dcu 2 , (32) 

where du> 2 is the standard line element on the unit sphere. The system is supplied with the 
boundary conditions 

lim A(t, r) = lim p(t, r) = X(t, 0) = , (33) 



3 We adopt the same notation as in the previous sections, although this differs from the standard notation for 
the Einstein-Vlasov system. 
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which define the asymptotically flat solutions with a regular center, and the initial condition 

< f(0,x,v) = r(x,v), f n (Ax,Av) = r(x,v) , VA G SO(3) . 

The reader is referred to [TH section 1.1] for a detailed derivation of the system. Throughout this 
paper we assume that / is a regular solution of ([30|) - ([3Tj) in the sense defined in [20]. In particular, 
f(t,x,v) is C 1 and has compact support in (x, v), for t G [0,T], and for any T > 0. For regular 
solutions, the metric coefficients are C 2 functions of their arguments. We emphasize that the 
existence and uniqueness of global regular solutions to the Cauchy problem for the system ([5D| - 
([3"T]) is open for general initial data. We also remark that the equation 

A r + fi r = Anre 2X {h + p rad ) , (34) 

follows by f3Ta |) -([31b | : by 434]) we have A r + \i r ^ and so, by (|33|) . 

0^ A + /i > /x(0, *) • (35) 



The ADM mass (or energy) H and the total rest mass M of a solution to the spherically symmetric 
Einstein- Vlasov system are defined by 

H= f f vT+|up / dv dx , M = [ [ e x fdvdx (36) 

and are constant for regular solutionsQ. Related to the ADM mass we have the quasi-local mass, 
defined by 

r r 

m(t, r)=4ir / s 2 h(t, s)ds=-(l- e~ 2A ) , (37) 
Jo 2 

where we used (|31al) . Thus Iim, — >0 o m(t, r) = H. 

For later convenience, we recall that the non-zero Christoffel symbols for the metric (|3"2"]l are given 

by 

1 00 - Mt i 1 0a - Mr— , 1 afc ~ e A t ~ "a - ' 



1 00 - e Mr , 1 06 - A * r 2 

\ XjjXa X C / c X\)X \ X a 

1 ab — A r 5 1 °b — 

r J r \ 

Note also that \g\ = e 2A+2M is the determinant of the metric. 

The stress-energy tensor T^ v for Vlasov matter in spherical symmetry is given by 

T 00 = e~ 2fl h , T 0a = e- x -^q — , (38a) 

r 

T ab = fi -2A rad XX + 1 tan U _ £1 (38b) 

H 2 \ r 2 / 

and satisfies the conservation lawd 

V^F" = . (39) 



4 The other two conserved quantities, the linear momentum Q and angular momentum L, considered in the 
previous sections are identically zero in the present context by spherical symmetry. 



5 The identities are a consequence of the Vlasov equation alone, see 1151 . 
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3.1 Virial identities for time dependent solutions 



To begin with we derive an integral identity for the spherically symmetric Einstein- Vlasov system 
as we did in Lemma [2] for the Nordstrom- Vlasov system, i.e., using the vector fields multipliers 
method. Actually, the identity in Lemma[3]below is valid not only for the Einstein- Vlasov system, 
but for all matter models in spherical symmetry. This is due to the fact that equation (|39[) . which 
is the starting point for deriving the integral identity, must be satisfied by all matter models for 
compatibility with the Einstein equations. 

Multiplying the conservation law ([3^)1 by a vector field £ M , integrating on a compact spacetime 
region Q with piecewise differentiable boundary dfl and applying the divergence theorem we obtain 
the integral identity 



J^rjn d<jg = 



OP. 



f T^V^dg 
Jo, 



(40) 



where rj^ is the normal covector related to the boundary, J M = — is the current 

associated to the vector field ^ and V M £„ = d^ u — r£„£ ff is the covariant derivative of the vector 
field. Moreover dg is the invariant volume clement on the spacetime and da g the invariant volume 
element induced on <9fL 

Lemma 3. Assume that (/i, <7,p rad ,p tan ) satisfy the compatibility conditio-^ (|39|) . where T^ v is 
the stress-energy tensor ([38]). In addition, we assume that h(t, •), q(t, ■), p rad (i, ■), p t&a {t, ■), have 
compact support. Given any smooth function x(i,r) in W lt J^° and any solution of (f3"Tj) define 



I(t) 



Then the following integral identity is verified: 



X q(t, r) dx . 



dl 



e M-A p rad 



dx 
dr 



dx 

at 



dx 



(41) 



Proof. In ([4"0f we use 



£o = 0, 

6 = x(t,r)- 



After a long but straightforward computation we obtain 



dx 
dr 



x 



,-2A 



hfi r + 2q\ t e 



— A — [i 



e - 2X p iad x r - p t 



We will choose f2 to be the coordinate image of a cylinder [0, T] x B(R). In this fashion, we have 
that 

l-T 



T^V^dg = 



x\<R 



dx 
dr 



z»- x X (hpi r +p™ d \ r 



dt 



dx dt. 



Now we compute the corresponding boundary integral in (j4H)l . First, the current reads 



qxe 



J a = e - 2A p rad 



X 

X(r) — 
r 



Next we write dfl — A\ U A2 U A3, where 



(42) 



3 In the case of a perfect fluid, the compatibility condition (139 I t is the system of Euler equations. 
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• Ai = {t = T, \x\ < R}. The outer unit normal is e^ r ' T ^dt; the induced metric is e 2A(r ' T) dr 2 + 
r 2 duj 2 , the volume element e 2A ( r ' T ) dx. 

• A2 = {t = 0, \x\ < R}. The outer unit normal is — e^ r '°^dt; the induced metric is e 2X ^ r, °^ dr 2 + 
r 2 duj 2 , the volume element e 2A ( r '°) dx. 

• A 3 = {0 < t < T, \x\ = R}. The outer unit normal has the form —e x ^ t,R ^^dxi. The metric 
is ds 2 = —e^^'^dt 2 + R 2 dw 2 , the volume element e 2 ^ t ' FCl dSndt, where dSn is the surface 
element on the sphere of radius R. 



Summing up we get 



da g = / q(T, r) X (T, r) dx - / q(0, r) X (0, r) dx 
an j\x\<r j\x\<r 

\ P tan (t, R) X (t, R) e ^,R)-Kt,K) dSR dt _ 

J\x\=R 



(43) 



Having assumed that the matter quantities are compactly supported in the variable r, the boundary 
integral vanishes in the limit R — > 00, whereas the other integrals remain boundeoj Thus in the 
limit we obtain 







/ qxdx 


T = fJ [ 


Jr 3 


JO Jr 3 . 



which is the integral version of (|4ip . 



dt 



2xA ( 



dx dt , 
□ 



We shall now derive two particular cases of the identity (|4"Tj) . First let us choose 

X = e 2X F(r) , 

for a smooth function F, We have d t x = 2X t e 2X F(r) and then dtX — ^X^t — 0. In this way equation 
(|4Tj) implies 













/ qe 2X Fdx 


T =fJ 


e M+A 


p Tad F' - 


f F 


Jr 3 


Jo Ji 


t 3 







„rad 



df . (44) 



Note now that, using flM), (|3TE|) and (J37J), 



A r p rad - Vr = P (Ar + Mr) - MP + A) = (A r + fX r ) P 



rad MrC 



-2\ 



Airr 



4 7rr 3 



(A r + /i r ). 



Then (|4"4"1) becomes 









/ qe 2X Fdx 


T = l T I e» +x 


p™ d F 


Jr 3 


Jo Jr 3 





,F _ F m(\ r + /y) 
r r 47rr 2 



dx dt . (45) 



In the integral in the right hand side we use that 



^+\ F (Ar + ^r)m 
ip" r Anr 2 



T f 00 de x+f " F 



dr 



-m dr dt 



7 Of course the compact support condition can be replaced by a suitable decay assumption. 
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T 

JO 



This leads to 



qe 2X Fdx 



Finally for F(r) = r we obtain 



±(^)e^drdt-H T (^^ 
dr \ r I \r— >oo r 



F(r) 



HT I lim 



r»H r-,1 tan ^ , ^ fXl d ( F 

p™ d F' +p tan - + /i— + - — ^— - 
r r 47rr z dr \ r 



J 



1 T 



qe 2X r dx 



- LI T + I / e x+f * (p rad + p tan + h) dx dt . 
Jo Jr 3 



dx dt 



(46) 



(47) 



3.2 Virial inequalities for steady states 

The existence of steady states solutions to the Einstein- Vlasov system is well understood, we refer 
to [3 [16] and the references therein. The identity (|47| restricted to steady states imply 



H 



e A+ "(p tan +p rad + h)dx 



(48) 



Remark 3. The fundamental identity (|48]) can be proved directly using the Einstein equations 
for static spherically symmetric spacetimes, see [3]. Our derivation has two advantages. Firstly, we 
obtained (|48|) as a special case of a more general identity which holds for time dependent solutions, 
see LemmaEl Secondly, the technique of the vector fields multipliers, which we used to derive (|48|) . 
can also be used on spacetimes which are not spherically symmetric and therefore our argument 
could be useful to prove generalizations of (|48[) for solutions with less symmetry. 



This identity leads naturally to a bound on the central redshift 

Z c = e -"<°> - 1 € [0, +oo) 

in terms of the mass-energy of the static solution. We consider only static solutions of the spheri- 
cally symmetric Einstein- Vlasov system. 

Proposition 1. Let f be a static solution of the spherically symmetric Einstein- Vlasov system 
with compact support. Then the following inequality holds true 



< 



ifH< M 
■2,,-u if H > M 



Z„ > 



H 



M 



(49) 



Proof. Since /x is increasing, fi(r) > fi(0) and so 



e x+fl (p Iad 



+ h) > e"<°> / e x h > Me^ . 



Using this in (|55]) gives 



H 

M 



(50) 
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Moreover 

JR 3 \ V 1 + M J ^R 3 V 



Thus, since A + [i > ^i(O) and < 1 < e A , 

e A+ f(p rad + p tan + h) dx > e^ 0) (2H - M) 



and so by (|48|) . 

The result follows from (|5D|) and (f5Tj) and taking into account that 

H H 

< 



M ~ 2H- M 

is satisfied in the case i < S < 1. □ 



3.3 Shells and Jeans' type steady states 



Let R be the radius support of the steady state. Using that [i is negative and increasing and that 
the steady state matches the Schwarzschild solution at r = R we obtain the bound 



e M0) < e „(«) = sjl-^-. (52) 

The inequality (f52"j) can be combined with (|4U)) to obtain an upper bound on e^^ ^ in terms of R, 
H and M. 



Now, we consider briefly an important class of steady states, namely the Jeans type steady states, 
see [28]. For these steady states the distribution function / has the form 

f(x, v) = tP(E, F) , where E = e^y/l + \v\ 2 , F = \x A v\ 2 . (53) 

Since the particles energy E and the angular momentum F are conserved quantities, the particle 
density (f5"3")) is automatically a solution of the Vlasov equation. The existence of Jeans type steady 
states is then obtained by replacing the ansatz / = ip(E,F) into the (time independent) Einstein 
equations and proving existence of global solutions for the resulting system of ODEs. We refer 
to [22] where this procedure is carried out for a large class of profiles i/j; moreover the Jeans type 
steady states constructed in [2 2) all have compact support and satisfy that 

3 E e (0, 1) such that ip = , for E > E Q . (54) 

Thus Eq is the maximum particle energy in the ensemble. The property (|54|) is necessary in order 
that the distribution function (|53[) be asymptotically flat and with finite energy. For Jeans type 
steady states one obtains a new estimate on e^ - 1 in a straightforward way: 

H=( { v /l + \v\ 2 fdvdx= ( ( e^^eV 1 + M 2 e A / dv dx < —LM , 
Jrz Jr3 Jr3 Jr3 e^> 



1G 



whence 

< E ^ . (55) 

In fact, for Jeans' type solutions we have [55] 



and thus, combining (|55[) with (|52[) wc conclude that for Jeans' type steady states the inequality 



e Mo)< e M*) =nim a,^l Ji 2/7 



H V i? 



hods. 



Consider now the case of a static shell. Let / be a static shell solution of the spherically symmetric 
Einstein- Vlasov system with inner radius R\ and outer radius Ri- Using (|31bj) we can write /i(0) 
as follows 



/u(0) = - [°° e 2A (^+4^' ad ) 



1 / ^ + 4^ rad )dr. 



o 1 — 2m /r \r 2 

By Buchdahl's inequality ([5]), the identity ([57)1 and the bound p rad < h, we obtain 



> 



/•"X 


/m 


■f Airrp lad ^ dr 


poo 

9 

J Ri 


(- 

\r 2 


+ Airrh^J dr 


9H 

i?7~ 


9 


f°° 2 

/ Airr h dr 
'Ri 


18H 






Ri ' 







Now, we use the upper estimates on /i(0) of the Proposition [T] to find 

/i(0) = In 

Combining both estimates we obtain 

18H 

Ri < 





1 < In | 


f 1 ^ 


(z c + ij 




M_i! 


+ 1 / 



In (If -1|+1) 



i.e., the inner radius of a static shell with given ADM energy and rest mass cannot be arbitrarily 
large. 
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